Path Crossing Exponents and the External Perimeter in 2D Percolation 
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2D Percolation path exponents xf describe probabilities for traversals of annuli by £ non-overlapping 
paths, each on either occupied or vacant clusters, with at least one of each type. We relate the 
probabilities rigorously to amplitudes of 0(N — 1) models whose exponents, believed to be exact, 
yield xf = (£ 2 — 1)/12. This extends to half-integers the Saleur-Duplantier exponents for k — 1/2 
clusters, yields the exact fractal dimension of the external cluster perimeter, Dep = 2 — x% — 4/3, 
and also explains the absence of narrow gate fjords, as originally found by Grossman and Aharony. 
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The fractal geometry of critical percolation clusters 
has been of interest both for intrinsic reasons and as 
a window on a range of phenomema. It is character- 
ized by fractal dimensions of various sets jj],U, e -S-' °f 
the connected clusters, their backbones, the sets of piv- 
otal (singly-connecting) bonds, the clusters' boundaries 
(hulls), and their external (accessible) perimeters. A set 
S is said here to be of fractal dimension Ds if the den- 
sity of points in S within a box of linear size R decays as 
R~ xs , with xs = d — D$ in d dimensions. 

For two-dimensional (2D) independent percolation, 
many of the fractal dimensions have been found exactly 
[^-^| , though most of these values have not yet been es- 
tablished at a rigorous level. In several cases, Saleur 
and Duplantier (SD) H identified the co-dimension x$ 
with the exponent x\~ which describes the decay law 
P<r w (r/R) x * for the probability (P£) that in an annular 
region D(r, R) the small circle of radius r is connected to 
the outer one, of radius R >> r, by k different clusters 
of occupied sites (or bonds). SD utilized the observa- 
tion that the statistics of the 2k boundary lines of the 
connected clusters correspond to those of loops in some 
well recognized models: the Q = 1 Potts model (at its 
critical point) for the bond percolation model and the 
0(N — 1) loop model of Domany et al. JtJ (at its low 
temperature phase) for site percolation on the triangular 
lattice. Using the "Coulomb gas" representation for the 



O(N) 



corresponding £-line exponents, 
for both models the values, expected to be universal, 



SD obtained §] 



Xu 



,°ir i} = - i)/i2 , 



(1) 



where k clusters correspond to £ — 2k lines in the loop 
model. 

Among the noteworthy applications of the above for- 
mula are the "hull dimension", i.e., the dimension of the 
cluster's perimeter, 



and the dimension of the set of "red" (singly connect- 
ing) bonds, which are pinching points between two large 
clusters: 



D sc = 2-x c 2 = 3/4 = l/v, 



(3) 
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where v is the correlation length exponent H, in agree- 
ment with previously derived values ^). However, some 
well known percolation dimensions have eluded this exact 
approach: the dimension Dep of the external (accessible) 
perimeter (EP) or frontier of a cluster, first studied by 
Grossman and Aharony (GA) B, and the backbone di- 
mension. The EP of a cluster is the accessible part of the 
hull, which excludes deep "fjords" which are connected 
to the cluster's complement through very narrow pas- 
sages (or "gates" ) . The dimension of the EP was found 
numerically to be Dep ~ 4/3 [jl(J. GA || also made 
the puzzling observation that, while typical clusters do 
show many fjords with only a narrow passage to the com- 
plement, once one fills in fjords with passages of width 
two or three lattice spacings - no fjords of broader mi- 
croscopic passages, and depth comparable with that of 
the cluster, are left. This is clearly visible in Fig. 6 of 
the second Ref. [|J . Both of these observations make the 
EP look very similar to self-avoiding walks (SAW's). Al- 
though there appeared conjectures attempting to make 
this relation quantitative j^J^] , the connection was never 
elucidated. 

In this Letter we report on a resolution of these issues 
through analysis of the path crossing probabilities. 
i. Basing the relation of percolation exponents with 
the 0(N = 1) exponents on a somewhat different footing 
than that used in Ref. M, we extend the list of exact 
values proposed for critical percolation in 2D. Instead of 
focusing on entire clusters, we consider the probability 
Pj '(r, R;ti, . . . , Tg) that the annulus D(r, R) is traversed 
by (at least) I non-overlapping connected paths, which 
are "monochromatic" in the sense that each consists of 
either occupied sites ("color" tj = +) or vacancies (tj = 
—). We rigorously prove O] that for color sequences 



1 



which include at least one of each type (±) the decay 
rates of the probabilities are color-independent, and are 
given by the 0(N = 1) exponents. Assuming the validity 
of the exact values for the latter, we find that 

PT(r,R;r 1 ,...,r e )^(r/R)^ (4) 
with the path crossing exponents xf satisfying 



■p _ 0(N=1) 

X X n 



(£ 2 - 1)/12. 



(5) 



Since the cluster exponents are x^ — xl^f., Eq. (||) may 
be viewed as an extension of the SD formula to odd val- 
ues of £, or half integer values of k, and to more general 
sequences r,- = ±. 




FIG. 1. A configuration with 1 = 2 paths traversing the 
semi-annular region D(r, R), one of occupied (checkered) and 
the other (dual) of vacant (gray) hexagons. The configuration 
has (£+1 = 3) 0(N = 1) boundary lines (drawn thickly); 
the states of the hexagons not adjacent to these lines are left 
unspecified. 



ii. Using the newly acquired values we explain some 
of the quantitative and qualitative features of the EP of 
critical clusters mentioned above: its dimension, which 
we identify as 



D 



EP 



= 2 



o(i) 



4/3 



(6) 



and the interesting fact that - unlike the hull - the EP 
appears to be self-avoiding on the macroscopic scale. 

in. We consider also the analogous boundary or 
"surface" exponents, which describe the probability 
Pj (r, R;t\, . . . ,Te) that, within the upper half space, a 
semi-annular region D(r, R) is traversed by £ paths (see 
Fig. 1). For the exponents defined by 



PT(r,R;n,..., n )*(r/Ry 



we find 



X p 



~0(N=1) 



(£ + l)£/6. 



(7) 



(8) 



In this case the relation is valid with no restriction on 
the color sequence t; however there is a shift: £ crossing 



paths correspond to (£ + 1) 0(N — 1) lines. Thus, with 
odd £, one recovers the cluster boundary exponents = 
x? k _ x = fc(2fc-l)/3, asinRefs. (l2g||. 

Before we turn to describe the arguments for the exact 
values of the path exponents, as provided by Eqs. (||) and 
(^), let us present their implications concerning the di- 
mension and shape of the external perimeter. Each point 
on the accessible EP is next to the end of three paths of 
lengths comparable with the diameter of the cluster - a 
path of occupied sites and in addition two distinct dual 
paths of vacancies (Fig. 2a), which guarantee that the 
point is not within a fjord of narrow opening (both paths 
must be able to exit the fjord via the narrow gate). This 



yields D EP = 2 — x 
the numerical results ll( 



4/3, in excellent agreement with 





a) 

FIG. 2. Paths, and dual (dotted) paths characterizing a) 
the external perimeter (I = 3) and b) a gate (£ = &). 



GA's observation concerning the fjords is particularly 
striking from the perspective of the scaling limit, for 
which one sends the lattice spacing to zero while keep- 
ing the sight on the curves observed on the macroscopic 
scale. (The limit can be constructed using the analy- 
sis of Ref . HH , which implies that the cluster hulls and 
EP's can still be described by means of Holder continu- 
ous random curves.) While the EP is self- avoiding on the 
lattice scale, like the hull it could have close encounters 
which appear as self-intersections when viewed from the 
macroscopic perspective. Yet such close encounters are 
not observed. Also this puzzle is explained by the gener- 
alized path statistics: the occurence in an L x L box of 
a cluster with a fjord of depth R — sL and neck width 
r = eL requires there being six paths, two pairs of triplets 
as used in the derivation of Eq. , which meet in a re- 
gion of size r and avoid each other up to a radius R (Fig. 
2b). The probability of finding such six paths scales as 
£ -d x (i^e = Q(e x e- d ) where d = 2. Equation §) 
yields the exponent value x§ = and hence the proba- 
bility for a randomly picked configuration to exhibit such 
a gate tends to zero, in the situation where s is fixed at 
some non-infinitesimal value < s < 1, and e — > 0. The 
crucial point here is that the fractal dimension of the set 
of these gates is negative: 



Dr 



< 



(9) 



2 



This explains the asymptotic absence of macroscopic size 
fjords of neck width r anywhere between few multiples 
of the lattice spacing and eL, for any fixed e < 1. The 
above argument also implies that the EP will not exhibit 
peninsulas with narrow isthmuses, even though that con- 
dition was not built into the construction. 

It is of interest to recall here the suggestion which was 
made on a theoretical as well as a numerical basis [|l5],|l6| , 
that percolation's hull and EP dimensions coincide with 
the dimensions of polymers, respectively at the 0-point 
(the onset of collapse), or in the SAW state: 



D 



H 



Do 



Dep — Dsaw 



(10) 



It is natural to conjecture that in the scaling limit the EP 
coincides, in its local statistics, with a SAW. This may 
appear to be in conflict with the a-symmetry between the 
two sides of the EP, however the absence of peninsulas in 
the scaling limit suggests that the symmetry is restored 
asymptotically. 

The above results (||-||) are based on a rigorous relation 
of path crossing probabilities with amplitudes of a loop 
model which we shall now define, and on known exact 
values for the exponents of the latter (which still remain 
to be proven at a rigorous level). The arguments, which 
will be presented more completely in ref. fPH , are for- 
mulated for the special model of independent site (i.e., 
hexagon) percolation on the triangular lattice. For rea- 
sons of universality one may expect the conclusions to 
apply to other 2D percolation models, e.g., the bond 
percolation, and also to the statistics of the connected 
clusters of (+) or (— ) spins in the 2D Ising model at all 
temperatures above T c . 

The loop-model configurations, T, are collections of 
nonoverlapping loops and lines, in suitable subsets of the 
plane, which are allowed to have end-points only within 
prescribed regions. The weight of a configuration is 



W(T) 



(11) 



with Afp the number of closed lines (or "polygons" ) and 
A/g their total length (the number of bonds). For the 
particular case of hexagon percolation discussed here the 
fugacities are K = 1, N = 1, i.e., W(T) = 1 for all T §. 

A probability distribution of loop and line configu- 
rations in a prescribed region is defined by means of 
the weights W(T)/Z, with Z a suitable normalizing fac- 
tor. Let now P^ N \r, R) denote the probability that 
such a system of lines with no end points in the annu- 
lar domain D(r, R) contains at least £ lines traversing 
D(r,R). For a representation of the surface exponents, 
we also let P^ N \r, R) denote the corresponding event 
with the lines restricted to lie in the upper half plane, 
taken here with the "free boundary conditions" . A close 
variant of the quantity P^ N ^ (r, R) is the O(N) ampli- 
tude G^ NS> (r, R) which is defined as the sum over £ lines 



{71, . . . , j£ } spanning the annulus D(r, R) of the probabil- 
ity w(7x, • • • , 7i S r) that the lines are included in T. For 
N = 1, that probability reduces to the local expression 
H0: 



w(7i 



(12) 



with Af H (71, . . . , 7^) the number of hexagons touched by 
the lines, £(71, ■ ■ ■ , 7^) the number of line clusters - two 
lines being regarded as in the same cluster if they touch 
a common hexagon, and # defined as taking the value 
if the £ lines leave room for another curve to traverse the 
annulus and 1 otherwise. The amplitudes then read 



G° {N=1 \r,R) = 



E 



w(7i) ■ ■ ■,7<), 



(13) 



the sum running over sets of I nonoverlapping lines which 
traverse D(r,R). It can be shown that the probabilities 
and amplitudes agree to the leading order pi|: 



p; 



0(N=1) 



(r,R) = 1 



o(-))G^ N=1 \r,R). ,1.4) 



"Coulomb gas" and Bethe Ansatz methods P,p|p7|-]l9|] 
yield the conclusion that the loop model amplitudes, 
and thus also the probabilities, decay by power laws, 
G° (1) (r,i?) « (r/ R yr\ and G° (1) (r,i?) « (r/R)*? W , 
with the exponents taking the values given in Eq. (|5|) 
and Eq. (||). Our results rest now on the fact that the 
0(N = 1) line probabilities are of the same order of mag- 
nitude as the path crossing probabilities. Their compat- 
ibility is expressed in the following statement. 

Proposition In the site percolation model on the trian- 
gular lattice: 

1) For any "color sequence" {tj — ±}j =1 which in- 
cludes at least one of each kind (+ and 



Pf{r,R;r u ...,T t ) 



>0(N= 



~ 1} (r,R), 



(15) 



where A < B means that there are constants < c\ , C2 < 
00 with which c%A < B < c-iA uniformly in r and R. 
2) The surface probabilities satisfy 



PJ'(r,R;T 1 , . . . ,n) 



= P%r i \r,R), 



(16) 



without any restriction on the color sequence r . 



Let us outline here the proof, whose details will be 
spelled in ref. |0. The simplest case of the above rela- 
tion is in the example of the half-disk amplitude with al- 
ternating color paths (as in Fig. 1), which corresponds to 
Pf{r, R; +, — , +,—,.. .). Equation ( |l6| ) holds there since 
the statistics of the boundary lines is given exactly by the 
0(N = 1) loop model. The result is then extended by 
establishing independence on the color sequence. This is 
done by successively conditioning on the suitable "right- 
most path" and flipping the site variables left of the line. 
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Thus use is made of the Markov property combined with 
the spin flip symmetry, which are enabled by the indepen- 
dence and the self duality of the site percolation model 
on the triangular lattice. The argument is a bit more 
involved in the case of the full disk. There we need to 
have at least one traversing boundary line, which we em- 
ploy to slit the annulus. The previous argument is then 
applied to the resulting simply connected domain. Equa- 
tion ( |l5| ) reflects the fact that the overcounting involved 
in the selection of the slit is by at most a finite factor. 

As noted in |20|pl| ], it is possible to obtain some se- 
lected path exponents by direct arguments. The values 
agree with the formulas given above. It is instructive to 
list specific values of xf = (£? — 1)/12: 



x% yields D H = Dg = 7/4. 
xf yields D E p = D S aw = 4/3. 
xl yields D sc = v~ x = 3/4. 
xf = 2 can be derived directly. 
xf > 2 implies that the EP is self-avoiding on 



1=2 

e = 3 
e = A 
i = h 
e = 6 

the large scale (Dg < 0). 

The relation (|5j> was not claimed for I = 1 , or for paths 
of a single color (xf). Concerning this let us note: 



It can be shown directly that x 1 



0(N=1) 



11 while 



x[ > |E2j. The path exponent is related to the cluster 
dimension; its value appears to be xf — 5/48 ||. 
- The case of £ — 2, with two paths of the same color, is of 
special interest since it relates to the backbone dimension. 
Numerically, xg = 0.3568 ± 0.0008 §|. 

For the surface exponents, xf = (1 + 1)1/6 of Eq. (||), 
we note that for £ — > oo xf /xf — ► 2, as it should; and 
1=1: xf = 1/3 is consistent with Cardy's equation 
for the crossing probability ||. 
£ = 2: x% = 1 can be derived directly. 



2 is also directly derivable. 



1 = 1 

The last one 21 is related to a slit-disc exponent which 
is attributed to J. van den Berg in Ref. pc| ]. 

Finally, we note that the SD formalism also yields pre- 
dictions for the hull dimensions of Fortuin-Kastclcyn ran- 
dom clusters, describing the Q-state Potts model. These 
were recently confirmed in numerical simulations by Hovi 



and Mandelbrot |24|. In contrast, the values found in 
that work for the external perimeters do not agree with 
the generalizations of the SD formulas to odd I. The 
results presented here were derived only for site percola- 
tion. It would be interesting to see generalizations. 

This paper is dedicated to the memory of Tal Gross- 
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Henri Poincare, the Institute for Advanced Studies (MA 
and BD), and of Tel Aviv University (MA). It was sup- 
ported in part by the NSF Grant PHY-9512729 (MA), a 
grant from the German Israeli Foundation (AA), and by 
a grant to the IAS from the NEC Research Institute. 



[1] 
[2] 

[3] 

[4] 

[5] 
[6] 

[7] 

[8] 
[9] 

[10] 



[11] 

[12] 
[13] 
[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 



[22] 
[23] 
[24] 



D. Stauffer and A. Aharony, Introduction to Percolation 
Theory (Taylor and Francis, London, 1994). 

H.E. Stanley, in Percolation Theory and Ergodic The- 
ory of Infinite Particle Systems, edited by H. Kesten 
(Springer-Verlag, IAM Volumes 8, 1987). 
M. den Nijs, J. Phys. A 12, 1857 (1979); Phys. Rev. B27, 
1674 (1983). 

B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982); J. Stat. 
Phys. 34, 731 (1984); in Phase Transitions and Criti- 
cal Phenomena, edited by C. Domb and J. L. Lebowitz, 
(Academic, London, 1987), Vol. 11. 

H. Saleur and B. Duplantier, Phys. Rev. Lett. 58, 2325 
(1987). 

J. L. Cardy, Nucl. Phys. B240 [FS12], 514 (1984); J. 
Phys. A 25, L201 (1992). 

E. Domany, D. Mukamel, B. Nienhuis and A. Schwim- 
mer, Nucl. Phys. B190 [FS3], 279 (1981). 

A. Coniglio, J. Phys. A 15, 3829 (1982). 

T. Grossman and A. Aharony, J. Phys. A19, L745 (1986); 
ibid. 20, L1193 (1987). 

GA found D E p = 1.37 ± 0.03 for E\ (second neighbor) 
gates on the square lattice, and 1.34±0.03 for E2 and E3 
(third and fourth neighbor gates) on the square and tri- 
angular lattices || . Similar values were measured numer- 
ically by P. Meakin and F. Family [Phys. Rev. A34, 2558 

(1986) ], and by M. Rosso [J. Phys. A 22, L131 (1989)]; 
and experimentally by A. Birovljev et al. [Phys. Rev. 
Lett. 67, 584 (1991)] for invasion percolation fronts and 
by L. Balazs [Phys. Rev. E54, 1183 (1996) for Al pitting. 
M. Aizenman, B. Duplantier and A. Aharony, in prepa- 
ration. 

B. Duplantier, Phys. Rep. 184, 229 (1989). 
J. L. Cardy, J. Phys. A 31, L105 (1998). 

M. Aizenman and A. Burchard, Duke Math. J. (to ap- 
pear) . 

A. Coniglio, N. Jan, I. Majid and H. E. Stanley, Phys. 
Rev. B35, 3617 (1987). 

B. Duplantier and H. Saleur, Phys. Rev. Lett. 59, 541 

(1987) . 

M. T. Batchelor and H. W. J. Blote, Phys. Rev. Lett. 61, 
138 (1988). 

B. Duplantier and H. Saleur, Phys. Rev. Lett. 57, 3179 
(1986). 

M. T. Batchelor and J. Suzuki, J. Phys. A 26, L729 
(1993). 

H. Kesten, Commun. Math. Phys. 109, 109 (1987). 
M. Aizenman, in Statphys 19, edited by H. Bailin (World 
Scientific, Singapore, 1996); and in Mathematics of Mul- 
tiscale Materials, edited by K.M. Golden et al, The IMA 
Volumes in Mathematics, v. 99 (Springer, 1998). 
H. Kesten, Pro b. Theor. Rel. Field s 73, 369 (1986). 



P. Grassberger, cond-mat/9808095 

J. -P. Hovi and B. B. Mandelbrot (unpublished). 



4 



